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pL^ ■ Abstract. High-resolution spectra for 24 SMC and Galactic B-type supergiants have been analysed to estimate the contribu- 

tions of both macroturbulence and rotation to the broadening of their metal lines. Two different methodologies are considered, 
viz. goodness-of-fit comparisons between observed and theoretical line profiles and identifying zeros in the Fourier transforms 
of the observed profiles. The advantages and limitations of the two methods are briefly discussed with the latter techniques being 
, adopted for estimating projected rotational velocities (v sin f) but the former being used to estimate macroturbulent velocities. 

■ The projected rotational velocity estimates range from approximately 20 to 60 km s" 1 , apart from one SMC supergiant, Sk 191, 

with a vsini~ 90km s _1 . Apart from Sk 191, the distribution of projected rotational velocities as a function of spectral type 
I/") , are similar in both our Galactic and SMC samples with larger values being found at earlier spectral types. There is marginal 

evidence for the projected rotational velocities in the SMC being higher than those in the Galactic targets but any differences 
are only of the order of 5-10 km s~ , whilst evolutionary models predict differences in this effective temperature range of typi- 
cally 20 to 70 km s~ . The combined sample is consistent with a linear variation of projected rotational velocity with effective 
temperature, which would imply rotational velocities for supergiants of 70 km s" 1 at an effective temperature of 28 000 K 
(approximately B0 spectral type) decreasing to 32 km s _1 at 12 000 K (B8 spectral type). For all targets, the macroturbulent 
broadening would appear to be consistent with a Gaussian distribution (although other distributions cannot be discounted) with 
Oh an ~ half-width varying from approximately 20 km s _1 at B8 to 60 km s _1 at B0 spectral types. 
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' 1 . Introduction with the predictions (see Fig. 7 of Trundle et al.). Additionally 

^3 . j—. in most previous analyses, it has been assumed that the line 

. . . , Evolutionary models, such as those by Heser & Lanser (2000 1, u , • ■ • tl ■ . .. 4 . r ... 

1 -L . j b b r> broadening was predominantly due to rotation and if this were 

Meynet & Maeder (2000 1 and Maeder & Meynet (2001 1, pro- . ,. , . , t . , 

J ■ ' 1 1 ; v not the case the discrepancy between observation and theory 

vide predictions of the evolution of both surface chemical com- ,,, , 

r would become larger, 
positions and rotational velocities of massive hot stars as a TT ^ . . . , . .. ... . TTT _ 

„ . „ . „ , „, Howarth et al. ( 1997 1 investigated the line widths in IUE 



in 
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function of time. One common theme of these models is that an 
initially high rotational velocity (of the order of 300 km s _1 ) is 



spectra for a large sample of early-type evolved stars and noted 

that there were no narrow-lined objects. They interpreted this 

required to produce significant mixing of nucleosynthetically ., ... t , , ^ , 

^ r b b j j as evlc [ ence that there was another mechanism (hereafter desig- 

processed material to the surface. Although rotational veloc- . , t , , , , ,. 

, . , . . ° , nated as macroturbulence ), which was contributing to the line 

lties decrease as their progenitor OB -type stars evolve away t , r> . i j^Twii i. i j 

f tU m . v> , ■ * u a broadening. More recently Ryans et al. ( 2002 1 have analysed 

high signal-to-noise (S/N) ratio (of the order of 500) spectra 
of Galactic supergiants in order to try and estimate the rela- 
tive contributions of rotation and macroturbulence. They found 

^ , , , that although the line widths implied projected rotational ve- 

giants - Hilher et al. 2003 ; O-type dwarf s - Bouret et al. 2003 , lU rcn , ,„„, -i . 

° , — Jr ___ — locities in the range of 50 to 100 km s , their best estimates 

B-type giants and supergiants - Lennon et al. 2003 Trundle et „ . orvi -i -.u u 

/ ^ f. , . , — ', — , were typically in the range of 10 to 30km s , with all upper 

al. 2004 1 appear to be systematically too low when compared , . , ^, _i . . . 

rr J r limits being less than 60 km s . These results were lower than, 

Send offprint requests to: P.L. Dufton, e-mail: for example, the calculations of Meynet & Maeder ( 20(50) for 

p . cluf ton@qub .ac.uk solar metallicity objects that predicted rotational velocities in 



from the main sequence, B-type supergiants showing processed 
material at their surface are predicted to retain significant rota- 
tional velocities. However measurement of projected rotational 
velocities for a range of SMC early-type stars (O-type super- 
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the range of 50 to 200km s _1 . Even allowing for projection 
effects, the observed values appeared to be systematically too 
low. 

Recently there have been two non-LTE analyses of high 
quality spectroscopy of SMC B-type supergiants. Trundle et al. 
( 2004 1 analysed the spectra of 7 supergiants (and a giant) us- 
ing the unified code FASTWIND (Santolaya-Rey et al. H997i to 
determine their photospheric and wind parameters. Dufton et 
al. J2005t used the static code TLUSTY (Hubeny 1988, Hubeny 
& Lanz 1995 1 to analyse 9 supergiants. Two of these targets 
were in common with those of Trundle et al. and a compari- 
son of the two analyses showed excellent agreement. Here we 
analyse the spectroscopic data from these two SMC investi- 
gations using the techniques discussed in Ryans et al. (|2002 1 
and Simon-Diaz & Herrero I2QQ51 to estimate projected rota- 
tional and macroturbulent velocities. These techniques utilise 
different methodologies and hence allow us to investigate their 
strengths and limitations. We conclude that the Fourier trans- 
form technique of Simon-Diaz et al. provides the better esti- 
mates for the projected rotational velocity and have hence used 
it to reanalyse the Galactic stellar data presented by Ryans et 
al. d2002l . 

Our principle aim is to obtain reliable estimates of the stel- 
lar projected rotational velocities while allowing for the con- 
tribution of other broadening mechanisms. In particular, we 
will investigate whether the low estimates found by Ryans et 
al. (120021 are confirmed and replicated in our SMC samples. 
One advantage of studying targets in the low metallicity en- 
vironment of the SMC is that the presence and the amount of 
contamination of the stellar surface by nucleosynthetically pro- 
cessed material is far easier to characterise. This in turn makes 
it easier to search for any correlation between rotation and the 
amount of such contamination. 



2. Observational data 

Our observational material consists of high resolution spec- 
troscopy for 13 SMC supergiants obtained at the Anglo- 
Australian Telescope (AAT) and the ESO Very Large Telescope 
(VLT) and for 11 Galactic supergiants obtained at the 
McDonald Observatory and the William Herschel Telescope 
(WHT). 

Trundle et al. ( 2004 1 obtained spectroscopy of 7 SMC tar- 
gets at the VLT with UVES instrument in November 2001. Of 
particular relevance here are the observations obtained in the 
blue arm of the spectrograph with a wavelength coverage of 
3900 to 5000A. Trundle et al. rebinned their data to a pixel 
size of 0.2A (yielding signal-to-noise ratios ranging from 120 
to 190). Here we have used the unbinned spectra which have a 
spectral resolution, R— 40000. Two targets, AzV 104 (B0.5Ia) 
and AzV 216 (Bl III) have not been included in the current 
investigation — the former as the quality of the observational 
data was insufficient to accurately estimate the projected rota- 
tional velocity and the latter as it is not a supergiant. The VLT 
data are summarized in Tableland further details of the obser- 
vations and the reduction techniques can be found in Trundle 
et al. (2004). 



Spectroscopy for 9 SMC B-type supergiants (two in com- 
mon with the VLT observations) was obtained at the Anglo- 
Australian Telescope in September 1996, with the UCLES in- 
strument. There was complete wavelength coverage form 3900 
to 4900 A with a pixel size of ~0.08 A, a spectral resolution, 
R^20 000 and with signal-to-noise ratios ranging from 90 to 
130. Note that although the signal-to-noise ratios of the AAT 
data is lower than that of the binned VLT data, the differences 
in pixel size lead to the data being of similar quality (see, for 
example, Fig. 3 of Dufton et al. 2005 1. The stars observed are 
again summarized in Table |2] and further details of the reduc- 
tion techniques can be found in Dufton et al. ( I2005i . 

Ryans et al. (2002i has presented spectroscopy for the 
Galactic supergiants with signal-to-noise ratios ranging from 
250 to 700 and spectral resolutions of R-160 000 (McDonald) 
and R-40000 (WHT). HD 36371 (B4 lab) has been excluded 
from this investigation due to the difficulty of estimating the 
projected rotational velocity with the other targets being listed 
in Table 121 Full details of the observations and their reduction 
can be found in Ryans et al. 

3. Analysis and Results 

3. 1. Analysis methodologies 

We have considered two methodologies to investigate the 
broadening of the absorption lines in our supergiant sample and 
in particular to distinguish between rotational broadening and 
that due to other mechanisms that we characterise as macro- 
turbulence. These utilise the goodness-of-fit between observed 
and theoretical profiles and the shape of the Fourier transforms 
of observed absorption lines and are discussed in detail by 
Ryans et al. (|2002 ) and Simon-Diaz & Herrero ( 2006 1 respec- 
tively. Below we summarize the two methodologies. 

The methodology of Ryans et al. is, in principle, straight- 
forward. Equivalent widths of relatively strong, unblended 
lines in each star are measured, and then theoretical profiles 
of the same strength are generated using a non-LTE code (e.g. 
TLUSTY - Hubeny 1988 1 with appropriate atmospheric param- 
eters. It is assumed that the broadening of the intrinsic absorp- 
tion line profile is due to rotation (v sin z) and another mech- 
anism that is characterised as macroturbulence and is initially 
assumed to have a Gaussian distribution of velocities (with a 
half-width at which the profile has fallen to - of its maximum 
value given by v m ). Theoretical profiles are then convolved 
with these functions and compared to the observational data. 
The sum of the squares of the differences are computed, giving 
a measure of the quality of the fit at each point in (v m , vsin/) 
space. The main difficulty with this goodness-of-fit (GOF) ap- 
proach is that profiles generated with different (v m , vsin/) 
combinations are quite similar especially when one broaden- 
ing mechanism dominates and we return to this point in Sect 

ED 

By contrast the Fourier transform (FT) methodology dis- 
cussed by Simon-Diaz & Herrero (2006 1 allows the projected 
rotational velocity to be derived independently of any other 
broadening mechanism which may affect the line profile. A 
line profile can be considered to consist of the convolution of 
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Table 1. Comparison of the estimates of the projected rota- 
tional velocities found using goodness-of-fit (GOF) and Fourier 
transform (FT) techniques. For both cases, the ranges of the 
estimates and their medians are listed. Also tabulated are the 
maximum possible values of v sin i taken directly from Ryans 
et al. together with the number (n) of lines considered in the 
current comparison. 



Star 


Goodness-of-fit 


Fourier Transform 


n 




Range 


Median 


Max. 


Range 


Median 




HD 2905 


1-22 


20 


<60 


51-59 


52 


4 


HD 13854 


11-15 


13 


<40 


31-34 


32 


3 


HD21291 


5-8 


7 


<20 


16-18 


17 


2 


HD 24398 


4-19 


9 


<40 


34-39 


36 


4 


HD 34085 


4-15 


10 


<30 


20-26 


25 


6 


HD 38771 


5-21 


9 


<50 


46-55 


50 


4 


HD41117 


3-24 


14 


<40 


40-42 


41 


2 


HD 193183 


2^17 


20 


<50 


40-50 


45 


11 


HD 204172 


3-71 


57 


<60 


57-70 


62 


8 


HD 206165 


3^4-1 


11 


<40 


25^U 


35 


7 


HD 208501 


9-11 


11 


<30 


25-30 


28 


4 



natural, instrumental, rotational and macroturbulent broaden- 
ing profiles. As first discussed by Carroll J1933al IT933b > the 
Fourier transform of the rotational profile has zeros whose po- 
sition in frequency space depend on the projected rotational ve- 
locity. The FT method for the determination of v sin i (c.f. Gray 
I1973i is then based on the fact that in Fourier space, convolu- 
tions transform into products. Hence it is possible to determine 
the projected rotational velocity of a star once the position of 
the first of those zeroes is identified in the Fourier transform of 
the line profile. Although this methodology has not been widely 
used in the study of OB-type stars (e.g. Ebbets 1979 ), the recent 
study by Simon-Diaz & Herrero ( 2006 1 illustrates the strength 
of this method for determining projected rotational velocities 
in early type stars (see also Simon-Diaz et al. 2005 1. 

3.2. Comparison of methodologies 

We have used the high quality spectroscopy reported in Ryans 
et al. ( 2002 1 to compare the results obtained by the two method- 
ologies. We note that the FT approach does not directly provide 
estimates of the macroturbulence parameter (v m ) and hence 
the comparison is restricted to the estimates of the projected 
rotational velocities. In Tabled we summarize the estimates 
from the two methodologies with the goodness-of-fit results be- 
ing taken directly from Ryans et al. Note that for some lines, 
it was not possible to identify the zeros in the Fourier trans- 
forms and these cases have been excluded from the compar- 
ison. Consequently the GOF results listed in Table ^ differ 
slightly from those reported by Ryans et al. The only signifi- 
cant change is for HD 204172 where the exclusion of the Mg II 
line at 4481 A and the Si III line at 4552 A leads to an increase 
in the median from the GOF estimates. 

Inspection of Table ^ implies that the estimates from the 
GOF technique are systematically lower than those from the 



FT technique with the mean of the differences being 22±10 
km s ~ 1 . However the FT estimates are consistent with the max- 
imum values deduced from the GOF analysis. We believe that 
these differences reflect the difficulty in identifying and quan- 
tifying the relatively small amount of rotational broadening in 
the presence of significant macroturbulent broadening. 

With the GOF methodology, the profiles become effectively 
degenerate for small values of v sin i with the contour maps in 
the (v m , vsinz') plane having broad flat ridges along the vsin/ 
axis (see, for example Fig. 3 of Ryans et al. 2002 1. Hence esti- 
mates of the projected rotational velocities become sensitive to 
the observational uncertainties and the assumptions adopted. 
For example, a Gaussian velocity distribution has been as- 
sumed for the macroturbulence and if the distribution differed 
from this, it could lead to systematic errors in the estimation 
of the projected rotational velocity. Additionally the sensitiv- 
ity of the GOF approach to uncertainties in the observational 
data (such as blending and continuum normalisation) may be 
the cause of the wider range of estimates found in a given star 
compared with the FT approach. 

For slowly rotating stars, the Fourier transform of the 
broadening function will have its first (and subsequent zeros) 
at high frequencies. Additionally if the observed stellar absorp- 
tion lines are significantly broadened by macroturbulence most 
of their power will occur at low frequencies in the Fourier do- 
main. Since the noise in the spectrum transforms as a 'white 
noise' in the Fourier space (Smith & Gray 1976), then the un- 
ambiguous identification of the zeros in the Fourier transform 
becomes difficult or impossible when the v sin i is low and the 
effect of macroturbulence is important (this is analogous to the 
difficulty in the GOF approach of estimating projected rota- 
tional velocities when the line profiles are dominated by macro- 
turbulence). Indeed even for the high quality Galactic super- 
giant spectra it was not possible to estimate the projected ro- 
tational velocity from the Fourier transform for some spectral 
features in a given target. Tests undertaken by Simon-Diaz & 
Herrero (2006 1 for cases with small projected rotational veloc- 
ities and significant macroturbulence indicate that the former 
can be overestimated by 5 to 10 km s -1 . Although this would 
reduce the discrepancy with the estimates deduced from the 
GOF methodology, it would not eliminate it. 

However the advantage of this methodology is that if the 
zeros in the Fourier Transform can be unambiguously identi- 
fied the corresponding v sin i estimates should be robust and in 
particular they are less dependent (compared with the GOF ap- 
proach) on additional assumptions or on uncertainties in the ob- 
servational data. Indeed indirect evidence for this is the smaller 
range in estimates found from different lines in a given star us- 
ing the FT methodology. 

3.3. Estimates of the projected rotational velocity 
and macroturbulence 

Given the discussion of the previous section, we concluded that 
the FT methodology provided the better estimates of the pro- 
jected rotational velocities and adopted this approach for both 
our Galactic and SMC samples. These estimates are summa- 
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Table 2. Median values for the macroturbulent (v m ) and 
projected rotational (vsiru) velocities deduced for the SMC 
datasets presented in Dufton et al. ( 2005 1 and by Trundle et al. 
(2004 1. Also listed are estimates for the Galactic supergiants 
discussed by Ryans et al. (|2002). The sources of the spectral 
types and effective temperatures are discussed in the text. 



Star 


Spectral 




Dataset 


v sin i 


Vm 




type 


K 




km s~ 




SMC 












AzV487 


BCOIa 


27000 


AAT 


50 


54 


AzV215 


BNOIa 


26750 


AAT/VLT 


66 


55 


AzV 242 


Blla 


21500 


AAT 


54 


43 


AzV78 


B1.5Ia+ 


21250 


AAT 


46 


36 


Skl91 


B1.5Ia 


21500 


VLT/AAT 


92 


49 


AzV 210 


B1.5Ia 


20500 


VLT 


37 


36 


AzV 462 


B1.5Ia 


19000 


AAT 


43 


26 


AzV 303 


B1.5Iab 


18000 


AAT 


42 


26 


AzV 472 


B2Ia 


19000 


AAT 


33 


27 


AzV 18 


B2Ia 


19000 


VLT 


33 


30 


AzV 374 


B2Ib 


18500 


AAT 


47 


20 


AzV 362 


B3Ia 


14000 


VLT 


38 


30 


AzV 22 


B5Ia 


14500 


VLT 


33 


23 


Galactic 












HD 204172 


B0.2 la 


28500 


WHT 


62 


62 


HD 38771 


B0.5 la 


27500 


McDonald 


50 


55 


HD 2905 


BC0.7 la 


24000 


McDonald 


52 


59 


HD 13854 


Bl lab 


23500 


McDonald 


32 


53 


HD 24398 


Bl lb 


23000 


McDonald 


36 


38 


HD 193183 


B1.5Ib 


22500 


WHT 


45 


40 


HD41117 


B2Ia 


19500 


McDonald 


41 


38 


HD 206165 


B2Ib 


20000 


WHT 


35 


34 


HD 34085 


B8Ia 


13000 


McDonald 


25 


28 


HD 208501 


B8Ib 


13000 


WHT 


28 


27 


HD 21291 


B9Ia 


11500 


McDonald 


17 


21 



rized in Table |2] together with spectral types and estimates of 
the stellar effective temperature. For the projected rotational 
velocities we have quoted medians as these will be less affected 
by any spurious measurements although very similar values 
would have been found if we had adopted the means of the es- 
timates. We also undertook an analysis of the SMC dataset us- 
ing the GOF methodology and will discuss this briefly in Sect. 
0] Spectral types and effective temperatures were taken from 
McErlean et al. {1999 1 for the Galactic sample and Trundle et 
al. (|2004i and Dufton et al. (I20051 for the SMC targets. Note 
that the three studies utilised different non-LTE codes and dif- 
ferent physical assumptions about e.g. the stellar winds and 
hence there may be systematic differences between the effec- 
tive temperature scales. For two SMC targets, effective tem- 
perature estimates are available from both Trundle et al. and 
Dufton et al. and they differ by 500K and 2000K (in Table 
we adopt the means of these values), which may be indicative 
of the magnitude of any systematic differences. 

The FT technique does not provide a direct estimate of the 
magnitude of other broadening mechanisms (that have here 



been characterised as macroturbulence) although this will af- 
fect the overall shape of the Fourier transform. Hence we have 
used the GOF technique to estimate this quantity. Our approach 
was for each star to adopt the projected rotational velocity esti- 
mate listed in Tableland then to re-analyse the relevant lines 
to determine the macroturbulence (v m ) that best fitted the pro- 
file. These estimates are also listed in Table [2] where we have 
again adopted the median of the individual estimates. 

The uncertainties in our estimates of the projected rota- 
tional velocity and macroturbulence will depend on several fac- 
tors including the quality of the observational data and the mag- 
nitude of the quantity being estimated — for example, an es- 
timate of a small macroturbulence in a target with significant 
rotational broadening will have a relatively high degree of un- 
certainty. In Table we list the range of v sin i estimates for 
each of our Galactic targets with a similar spread being found 
for our SMC targets. These translate into a typical uncertainty 
in the median value of the v sin i estimates of 10%. The spread 
in the estimates of the macroturbulent velocity for a given tar- 
get is similar and implies a typical uncertainty in the medians 
of 10%, apart from the smallest values (v m ^ 30km s _1 ) where 
the uncertainty may be larger. 

4. Discussion 

4.1. Macroturbulence 

The macroturbulence has initially been assumed to be isotropic 
and to follow a Gaussian distribution. The good agreement 
between the theoretical and observed profiles (see Fig. 1 of 
Ryans et al 2002 1 provides some evidence in support of this as- 
sumption. However we cannot rule out other distributions and 
in particular those that are qualitatively similar to a Gaussian 
distribution. As an example we consider anisotropic macro- 
turbulence characterised by a radial-tangential model (Gray, 
119751 11992> . In this model the velocity field is modelled by 
assuming that a fraction (A r ) of the material is moving radi- 
ally (but with a Gaussian distribution of velocities) with the 
remainder moving tangentially. The methodology was devel- 
oped to model the flow of convective cells observed in the Sun. 
Hence it may not be appropriate to the targets considered here 
and should be considered as illustrative. As an example we 
have fitted the Sim line at 4552 A in HD 206 165 adopting a 
radial-tangential model with A r =l.Q and 9 r =55 km s _1 (i.e. 
assuming a complete radial flow). However we note that as- 
suming that both the radial and tangential velocities are equal 
{A r = A t , O r =Ot=55 km s" 1 ) would have yielded an effec- 
tively identical profile. In addition, we have also calculated a 
profile for an isotropic macroturbulence with w m =35.5 km s _1 . 
The results are shown in Fig.^and both approaches are in good 
agreement with the observational data, but the radial tangential 
profile requires a larger macroturbulent velocity to be invoked 
(by approximately 20 km s _1 ). Hence we conclude that it is not 
possible with the current dataset to investigate the exact form 
or degree of the anisotropy of the macroturbulence. However 
what is clear is that additional broadening is present, which 
is consistent with a velocity distribution that is approximately 
Gaussian. 
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The finite spectral resolution of our observational data 
should not affect our estimates of the projected rotational ve- 
locity but could led to overestimates for the macro turbulence. 

60 I 1 1 1 r- 

Our lowest resolution spectroscopy (R^20 000) would have a 
- width, v; ~6 km s _1 assuming that the instrumental profile s 
was a Gaussian. Then if the macroturbulent velocity field also 1 60 
followed a Gaussian distribution, the estimates listed in Table|2] f 50 
would represent the instrumental profile and the actual macro- I m 
turbulent velocity added in quadrature. As an example, correct- 1 30 « 
ing our lowest estimate of v m ~ 21km s _1 for such an instru- f _ . • * 
mental profile would lead to a decrease of less than 1km s -1 . 3 

10 - 

For larger macroturbulent velocity estimates or better spectral 
resolutions, the corrections will be smaller and hence we have % i.e 14 

not attempted to apply them. 



1.8 2.0 2.2 2.4 

Effective temperature (K) / 10* 




Fig. 1. Observed and theoretical profiles for the Si III line at 
4552 A in HD 206165. The former is represented by a his- 
togram whilst the theoretical profiles assume an isotropic 
Gaussian velocity (dashed line) and a radial-tangential model 
with A r = 1.0 (dotted line). Note that both theoretical ap- 
proaches lead to a good fit with observation particularly in the 
core of the line, however the radially tangential profile requires 
a larger macroturblent velocity. 



Fig. 2. Estimates of macroturbulent (upper figure) and pro- 
jected rotational velocities (lower figure) plotted against effec- 
tive temperature. The Galactic supergiants are represented by 
squares and the SMC targets by circles. For the projected rota- 
tional velocities, the linear least squares fit (excluding Sk 191) 
is shown as a solid line, whilst the dashed line would repre- 
sent the rotational velocity assuming a random distribution of 
projection angles. 



Inspection of results summarized in Table [2] implies that 
the amount of macroturbulent broadening in the Galactic and 
SMC supergiants are comparable. To illustrate this in Fig. |2] 
our best estimates of the macroturbulent velocities, v m , are 
plotted against the effective temperature for both sets of su- 
pergiants. As discussed in Sect. 13.31 the latter have been taken 
from different sources and although there may be systematic 
differences between these scales, we believe that this should 
not invalidate the comparison. 

In general, both samples show a similar trend of macrotur- 
bulence with effective temperature and would imply that this 
is an intrinsic property of the star. Additionally as discussed 
previously the effect of this broadening is significant and inval- 
idates the use of line widths as a measure of rotational broaden- 
ing for B-type supergiants. This is consistent with the analysis 
of IUE spectra by Howarth et al. J1997i . who found signif- 
icant spectral broadening in all their supergiant spectra. This 



they interpreted as 'confirming that an important line broaden- 
ing mechanism in addition to rotation must be present in these 
objects'. 

For one SMC target, Sk 191, a relatively large projected ro- 
tational velocity is found from the FT methodology and this 
was confirmed by an independent GOF analysis. As can be seen 
from Fig. [2] our estimate of its macroturbulent velocity is com- 
patible with the other targets having similar effective temper- 
atures. By contrast adopting a small rotational velocity would 
have led to an anomalously large estimate for the macroturbu- 
lence of approximately 80km s _1 . 

If the additional spectra line broadening is due to macro- 
turbulence, the estimated velocities, v m , are large and imply a 
highly dynamic atmosphere. Additionally B-type supergiants 
are also found to require significant microturbulent velocities 
that need to be incorporated into the profiles of the line absorp- 
tion coefficient, with typically values being between 10 and 20 
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km s _1 (see, Crowther et al. 120051 and references therein). If 
these represent small scale turbulence, then they and the values 
of macroturbulence found here imply that the photospheres of 
B-type supergiants are highly turbulent on all distance scales. 
In turn the adoption of a stable photosphere in both static and 
unified stellar atmosphere codes must be considered problem- 
atic although excellent agreement between theory and observa- 
tion has been found using such methods (Crowther et al. 2005 
Trundle et al. l2004l Dufton et al. 120051 . 

Assuming that our macroturbulent velocities are due to 
a large scale velocity field, we can estimate its kinetic en- 
ergy. As an example we consider three relatively hot targets, 
HD 38771, AzV215 and AzV487 with atmospheric parame- 
ters, T c ff~27000K and log g~2.90. The evolutionary models 
of Meynet & Maeder (2000 2001 1 then imply an evolutionary 
mass of approximately 40M Q . Additionally, the correspond- 
ing TLUSTY model atmosphere implies that the mass of a col- 
umn of unit area down to a Rosseland mean opacity of unity is 
10 kgm m~ 2 . Adopting this as the extent of the atmosphere 
sampled by the absorption lines implies a total atmospheric 
mass of 6 x 10 19 kgm and a typical macroturbulent velocity 
of 60km s -1 would then lead to a kinetic energy in the veloc- 
ity field of approximately 10 29 J. It should be noted that this 
is only an order of magnitude estimate as it incorporates sev- 
eral assumptions and simplifications. However the bolometric 
luminosity of these supergiants is approximately 6 x 10 5 L Q or 
2 x 10 32 J s" 1 and the kinetic energy of the macroturbulence 
would be therefore negligible when considering the energy bal- 
ance in the atmosphere. McErlean et al. ( 1998 ) discuss the ther- 
mal velocities and sound speed for a model with atmospheric 
parameters (T c g=27 500 K, log g= 3.0 and a microturbulence 
of 12 km s ) similar to those of the three targets consid- 
ered here. For hydrogen the thermal velocity is typically 20 
km s _1 and the sound speed is 15-20 km s . Hence the en- 
ergy of a macroturbulent velocity field would be greater than 
the thermal energy of the plasma, whilst its velocity would be 
supersonic. Assuming that this additional broadening is indeed 
due to a large scale velocity field, its origin and nature remain 
unclear. 



4.2. Projected rotational velocities 

The estimates of the projected rotational velocities, summa- 
rized in Table|2 imply that the values for both the Galactic and 
SMC samples are similar and this is supported by Fig.[2]where 
the estimates are plotted against effective temperature. If all our 
targets had the same initial mass and rotational velocity, the 
effective temperature would represent their age and evolution- 
ary status. From a comparison with the evolutionary models of 
Schaller et al. ( 1992) for the Galactic targets (see McErlean et 
al.. [T999l for further details) and Maeder and Meynet ( I200H for 
the SMC targets (see Trundle et al., 2004 1, we have estimated 
evolutionary masses. These fall in the range 20 to 4OM , apart 
from the SMC target, AzV 78, which has an estimated initial 
mass of approximately 50M Q . This significant range in masses 
coupled with the initial rotational velocity of our targets being 



effectively unknown implies that this figure should be consid- 
ered as primarily an empirical representation. 

Theoretical predictions imply that the SMC supergiants 
should have larger rotational velocities (see Fig. [3] and discus- 
sion in Sect. l4.3t and there is some marginal evidence for this in 
Fig-El To investigate this we have independently estimated lin- 
ear least squares fits for the two samples (excluding the anoma- 
lous SMC supergiant Skl91 which is discussed below). The 
two fits are similar, differing by 4 km s _1 at T e g= 12 000K and 
9 km s _1 at 28 000K. Given the small size of the samples and 
the possibilities of systematic differences in the effective tem- 
perature scales of the two samples these differences cannot be 
considered significant. Hence we have also combined the two 
samples and the corresponding linear fit is shown in Fig. [2] as 
a solid line. The estimates of the projected rotational veloc- 
ity will be affected by the unknown angle of inclination (i). 
Assuming that this was randomly orientated would lead to an 
average value of sinj of j and the dashed line in the figure 
represents the least squares fit increased by the reciprocal of 
this factor. Although the sample is small it is encouraging that 
this appears to trace the upper envelope of of our estimates and 
would be consistent with these stars having an angle of inclina- 
tion close to 90°. Hence assuming that all B-type supergiants 
undergo the same evolution of rotational velocity (v ro t) with 
effective temperature our best estimate for that relationship is 
given by: 

v Iot = 3.4 + 2.4 x 10- 3 T cS (1) 

It is emphasised that this relationship should be treated with 
considerable caution for several reasons. Firstly our sample is 
relatively small with only 24 targets. Additionally it is unre- 
alistic to assume that every B-type supergiant will follow the 
same evolutionary history as this will depend on, for example, 
whether it is evolving on a blue loop having been a red super- 
giant and whether it is or was part of a binary system. Indeed 
the relatively tight relationship between projected rotational ve- 
locity and effective temperature (excluding Sk 191) in Fig.|2]is 
surprising. This is because a spread of estimates would be ex- 
pected due to both the random inclination of rotational axes 
and to variations in the initial rotational velocity. For example, 
targets with low v sin i estimates would be expected covering 
the full range of effective temperatures as these would have 
evolved from main sequence O-type stars with low rotational 
velocities. 

An explanation may lie in the biases in our observation 
samples. For example, two supergiants were excluded from the 
analysis as it was not possible to identify zeros in their Fourier 
Transforms, which could be due to them having intrinsically 
low projected rotational velocities. Hence our results could be 
biased to larger estimates, although as our initial sample size 
was twenty six targets any bias would be relatively small. 

For each part of our sample, the main criteria for target se- 
lection was as follows. For the Galactic supergiants, there was 
a large choice of possible targets and the final selection aimed 
to cover a range of spectral type for both la and lb luminosity 
classes. As such there is no clear source of bias in our selection. 
The VLT spectroscopy was obtained principally to investigate 
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properties of the stellar winds, whilst the AAT targets were cho- 
sen to cover a range of photospheric N/C abundances ratios 
(based on the analysis of McErlean et al. 1999) but with a bias 
to brighter targets for ease of observation. Hence both SMC 
samples may be biased to the highest luminosity objects and 
indeed only contain two objects with luminosity class lab or 
lb. Such samples might be expected to have intrinsically small 
rotational velocities due to their strong winds. However this 
would not explain either the lack of low vsin/ objects at the 
earliest spectra-types or the small scatter in the estimates. 

However given all the caveats, discussed above the relation- 
ship in equation^has some merit in indicating the magnitude 
of the rotational velocities that should be found in evolutionary 
models of blue supergiants. 

4.3. Comparison with evolutionary predictions 

Previously, several authors have estimated projected rotational 
velocities and compared these with the predictions of stellar 
evolutionary calculations. For the A-type supergiants (Venn 
I1999> . there appears to be reasonable agreement. However for 
O-type dwarfs (Bouret et al. 2003 1, O-type supergiants (Hiller 
et al. 120031 Crowther et al. 12002b . B-type giants (Lennon et 
al. 2003 1 and B-type supegiants (Trundle et al. 2004}, the ob- 
served values of v sin i appear to be smaller than those implied 
by stellar evolutionary calculations. However this comparison 
is complicated by comparing observed projected rotational ve- 
locities with predicted rotational velocities. Additionally it has 
normally been assumed that rotation was the dominant mecha- 
nism in the line broadening (with for example macroturbulence 
not being included) leading in effect to upper limits for the pro- 
jected rotational velocities. The advantage of the current esti- 
mates are that we have attempted to distinguish between the 
different broadening mechanisms and thereby have been able 
to estimate the actual values of projected rotational velocities. 
In turn, this has allowed us to deduce a relationship between 
the mean rotational velocity and effective temperature. 

In Fig. |3 we compare our results with those predicted 
by evolutionary models appropriate to both Galactic (Meynet 
& Maeder I20U31 and SMC metallicities (Meynet & Maeder 
I2004> . All the models illustrated have initial rotational veloci- 
ties of 300 km s _1 and have initial masses in the range of 20 to 
40 M Q . As the initial rotational velocities of our targets are un- 
known (and indeed likely to vary from target to target), this 
comparison should be taken as illustrative. However all our 
SMC targets show nuclearsynthetically processed material in 
their atmospheres and if rotation was indeed the main mecha- 
nism for this enrichment, it would appear that their initial ro- 
tational velocities must have been relatively large. Also shown 
are the measured projected rotational velocities and our esti- 
mate of the mean rotational velocity as a function of effective 
temperature taken from equation (1). For the Galactic super- 
giants, our results are in reasonable agreement with the pre- 
dictions for the 40M Q model but may be slightly lower than 
those for the 20M Q . However given the uncertainty in the inti- 
tial projected rotational velocities and the relatively small size 
of our sample, the agreement must be considered encouraging. 




LogPW) 

Fig. 3. Predictions of the evolution of the rotational veloc- 
ity of two models with metallicities appropriate to our Galaxy 
(Meynet & Maeder 2003 , initial masses: 20 M - narrow solid 
line, 40 M©- narrow dashed line) and two models representa- 
tive of the SMC (Meynet & Maeder l2004l 30 M Q - thick solid 
line, 40 M©- thick dashed line). All models had an initial ro- 
tational velocity of 300km s _1 . Also shown are the estimated 
projected rotational velocities for our Galactic (filled squares) 
and SMC (filled circles) samples. The SMC target Skl91 is 
represented by a filled circle enclosed in a diamond, for clarity. 
The dotted line represents our estimated mean rotational veloc- 
ity for both samples, excluding Sk 191 and corrected assuming 
a random inclination angle. 



The evolutionary models appropriate to our SMC sample pre- 
dict larger rotational velocities - for example at log T e g= 4.4 
dex, the projected rotational velocities for the 40M Q models are 
approximately 150 and 70 km s _1 at SMC and Galactic metal- 
licities respectively. Although the projected rotational velocity 
of Sk 191 is consistent with these predictions, those for most 
of the other SMC targets appear to be too small. This is con- 
firmed by the discrepancy between the theoretical predictions 
and our estimate of the rotational velocity at early B-spectral 
types (4.45 5s log T e ffS5 4.25), although the agreement may 
be better at later spectral types. We note that tests undertaken 
by Simon-Diaz & Herrero (2006 1 indicate that for relatively 
large macroturbulences and small projected rotational veloci- 
ties, the FT methodology may overestimate the latter. However 
allowing for this effect would increase the discrepancy between 
theory and observation for the SMC dataset. 

A comparison is also possible with the models of Heger 
& Langer (2000 1. Unfortunately they do not explicitly tabu- 
late rotational velocities as their models evolve from the main 
sequence. However their calculations imply that massive stars 
should evolve to higher effective temperature after their red 
supergiant stage (the so called 'blue loops'). Heger & Langer 
predict that such objects should have typical rotational veloci- 
ties of 50km s _1 and that their maximum effective temperature 
should be approximately 15 000K. For the late B-type objects, 
our estimated rotational velocities are smaller but compatible 
with these predictions. 
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Projected rotational velocity (km/s) 

Fig. 4. Projected rotational velocities plotted against nitrogen 
abundances for the SMC targets. The cross (+) is for the tar- 
get, AzV487, for which only an upper limit for the nitrogen 
abundance could be estimated. 



The theoretical predictions for objects evolving from the 
main sequence also imply that the amount of processed ma- 
terial mixed to the surface should depend on the initial stellar 
rotational velocity. For example Maeder & Mey net (2001 1 con- 
sidered 20M Q evolutionary models with an SMC metallicity 
and initial rotational velocities ranging from to 400km s _1 . 
At the end of the hydrogen burning phase, these models showed 
changes in their N/C abundance ratio, which were negligible 
for the non-rotating model but increased to a factor of approxi- 
mately six for the fastest rotating model. More recently Meynet 
& Maeder (2004i have discussed models with different initial 
metallicities. For models with significant rotational velocities, 
they find increases of approximately a factor of ten in the ni- 
trogen abundances as the stars evolve into blue supergiants. 
Additionally Howarth & Smith (2001 1 have discussed the pro- 
jected rotational velocities for ON and O-type main sequence 
stars. The cumulative probability function of the former im- 
plied that they were rotating more quickly with a Kolmogorov- 
Smirnov test indicating that the differences were statistically 
significant. Hence for these O-type stars (that will evolve into 
B-type supergiants), such a correlation between nitrogen abun- 
dances and rotational velocities appears to be present. 

We have searched for such a correlation in our SMC sample 
(we do not have reliable nitrogen abundances for our Galactic 
targets) by plotting projected rotational velocities against nitro- 
gen abundances in Fig. |4] The latter have been taken directly 
from Dufton et al. (2005) and Trundle et al. (2004 1, although 
as discussed in Dufton et al. (2005 1, for two stars in common 
there was a difference of approximately 0.15 dex in the nitro- 
gen abundance estimates. Tests showed that allowing for this 
offset made no substantial difference and hence it has not been 
included in the figure. No positive correlation of nitrogen abun- 
dance with projected rotational velocity is observed. Indeed 
there is a suggestion that our most nitrogen enriched targets 
have the lowest values of v sin However our sample size is 
relatively small and hence we cannot preclude the possibilty 



that a correlation has been masked by different angles of incli- 
nation and differences in spectral type. 

Although our targets normally have relatively small rota- 
tional velocities, there is one exception in the SMC supergiant 
Sk 191. We believe that the result for this star is robust as it 
is based on two different spectroscopic datasets and is sup- 
ported by analyses using both the GOF and FT methodologies. 
Additionally whilst the metal line profiles of the other super- 
giants appear to be approximately Gaussian, those of Skl91 
have a profile closer to that characterising rotation (see, for ex- 
ample, Gray 1992 1. However the estimated atmospheric param- 
eters and chemical composition of Skl91 are consistent with 
those of the other targets in our SMC sample and hence we are 
unable to explain the cause of this enhanced rotation. 

Hence we would appear to have a situation where effec- 
tively all our SMC supergiants have processed material at their 
surface with most of them having relatively small rotational ve- 
locities. Recently, consideration has been given to the effects of 
magnetic fields (Maeder & Meynet 120031 120041 120051 . These 
would appear to suppress mixing of nucleosynthetically pro- 
cessed material to the stellar surface. However when magnetic 
and thermal instabilities are also considered, these favour the 
chemical transport of elements but also lead to larger rotational 
velocities during the main sequence stage of, for example, a 
15M Q model considered by Maeder & Meynet (2005 1. Hence 
although magnetic fields may influence the presence of nucle- 
arsynthetically processed material at the surface of B-type su- 
pergiants, current models do not appear to explain the relatively 
low rotational velocity estimates found here. 

5. Conclusions 

Our principle conclusions are as follows: 

1 . We have compared two methodologies for investigating the 
additional line broadening in spectra of B-type supergiants. 
We find that Fourier Transforms of the absorption line pro- 
files yield the more reliable estimates of the projected rota- 
tional velocity, whilst a 'goodness-of-fit' methodology can 
be used to estimate macroturbulence. 

2. We have analysed AAT/UCLES and VLT/UVES spec- 
troscopy of 13 SMC supergiants to estimate the relative 
contributions of macroturbulence and rotation to the line 
broadening. These have been supplemented with a re- 
analysis of spectra of 1 1 Galactic supergiants. 

3. The estimates of the macroturbulence in the SMC and 
Galactic supergiants show no significant differences. 
Additionally both datasets show a similar behaviour of the 
macroturbulence increasing with effective temperature. 

4. With the exception of Skl91, our estimates of the pro- 
jected rotational velocities show no significant differences 
between the Galactic and SMC targets, although the latter 
may be systematically larger by 5-10km s _1 . As for the 
macroturbulence, the estimates of projected rotational ve- 
locity increase with effective temperature. 

5. Assuming a random orientation of the axis of rotation our 
results are consistent with a linear behaviour between rota- 
tional velocity and effective temperature with value of v TOt 
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of 60 and 30 km s _1 at effective temperatures of 28 000K 
and 12 000 respectively. 

6. For our SMC sample, we find no positive correlation be- 
tween nitrogen abundances and projected rotational veloc- 
ity as might be expected from stellar evolutionary calcula- 
tions. 

7. Skl91 has a macroturbulence consistent with its spectral 
type but an anomalously large projected rotational velocity 
of 9210ns- 1 . 
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